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Abstract 


The fringe-type laser anemometer is analyzed 
using the Cramer-Rao bound for the variance of the 
estimate of the Doppler frequency as a figure-of- 
merit. Mie scattering theory is used to calculate 
the Doppler signal wherein both the amplitude and 
phase of the scattered light are taken into 
account. The noise from wall scatter is calculated 
using the wall bidirectional reflectivity and the 
irradiance of the incident beams. A procedure is 
described to determine the optimum aperture mask 
for the probe volume located a given distance from 
a wall. The expected performance of counter-type 
processors is also discussed in relation to the 
Cramer-Rao bound. Numerical examples are presented 
for a coaxial backscatter anemometer. 

Nomenclature 


a expected peak pedestal counts per 

Doppler cycle 

A expected peak pedestal counts per second 

c velocity of light 

<> expected value of 

JE electric field vectors 

E q electric field amplitude 

fg frequency difference between incident 

beams, Hz 

fp Doppler frequency, Hz 

F lens focal length 

h Planck's constant 

I(x,y) irradiance at point (x,y) 

k wave number ( = UD 

2 wave vectors of incident beams 
k $ wave vector of scattered light 

L r radiance of laser light reflected 

by wal 1 

L^ radiance of background light 

£ measurement set vector 

n. number counts in i th time interval 

n Q expected noise, counts per Doppler cycle 

n(x,y) noise defined by eq. (10) 

n k noise counts due to flux through k th 

area element (eq. (21)) 
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Nf number fringes between 1/e irradiance 

points 

N q expected value of noise from wall scatter 

and background 
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unit vectors normal to incident beams 
and in scattering planes 
pedestal signal defined by eq. (11) 
likelihood function 
integral of p(x,y) over aperture 
laser power 

unit vectors normal to scattered beam 
and in scattering planes 
signal defined by eq. (11) 
integral of q(x,y) over aperture 
signal counts due to flux through k ^ 
area element (eq. (21)) 
quality factor defined by eq. (22) 
distance from scatterer to observation 
point 

unit vectors normal to scattering planes 
signal defined by eq. (11) 
integral of s(x,y) over aperture 
signal counts due to flux through k ^ 
area element (eq. (21)) 

Mie scattering amplitude function 
(a = "s\ "p") 
signal-to-noise ratio 
time, sec 

time at peak of Doppler burst, sec 
visibility of Doppler burst 
variance of ( ) 

beam radius at 1/e irradiance point 

beam waist radius 

intrinsic impedance of free space 


t h 

a.j v unknown parameter 

3 half of beam crossing angle 

£ Fisher information matrix 

6 a (e) phase of scattering light (a = "s", "p") 

6 aB phase difference of light scattered from 

two beams (eq. ( B13) ) 
aA area element in aperture plane 

Af system noise bandwidth 

At time interval 

n quantum efficiency 

0 angular coordinate 
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scattering angles 
wavelength of laser light 
bidirectional reflectivity 
scattering cross sections defined by 
eq. (B13) 

standard deviation of estimate of 
Doppler frequency 
angular coordinate 

flow angle relative to fringe normals 
phase of Doppler signal 
solid angle subtended by pinhole image at 
point in aperture planel 


Introduction 


The laser anemometer (LA) has become a valu- 
able tool in turbine engine research, providing 
data that would be almost impossible to gather 
using conventional instrumentation. However, this 
application of LA has proven to be one of its more 
difficult. Measurements in turbomachinery com- 
ponents are typified by small passages and highly 
accelerated, high velocity flows. This leads to 
the requirement for small seed particles that will 
faithfully follow the flow. Unfortunately, small 
particles are weak light scatterers, which result 
in low signal levels. When coupled with the high 
level of noise caused by laser light scattered 
from surfaces near the measurement point, the 
signal-to-noise ratio (SNR) is often less than 
needed for accurate determination of the particle 
velocity. 


Limitations of the conventional dual-beam 
fringe-type LA have forced researchers to look for 
new techniques and new optical designs. 

In particular, the dual-focus (two-spot) LA 
was developed^/ for use in engine research. 

This system is more sensitive to small particles 
and has better stray light rejection properties 
than the fringe LA, However, the dual-focus LA 
generally is more expensive and requires more 
sophisticated signal processing than the fringe 
LA. Also, in highly turbulent flows, the data 
rate tends to be lower than the data rate achieved 
with the fringe LA. (This low data problem may be 
solved by use of elliptical spot systems such as 
described in ref. 2.) 


One technique that has permitted the frinqe 
UA to be successfully used in compressor and tur- 

tu use of f ^orescent seed par- 
ticles ) The laser-induced fluorescence from 
the particles is separated from the laser light 
scattered from walls by an optical filter. This 
allows measurements to be made close to surfaces 
Unfortunately, the fluorescent efficiency of the" 
dyes used decreases rapidly as the temperature of 
the flow increases above room temperature. Thus 
unless a high efficiency, high temperature fluo-’ 
rescent seed material is found, it appears that 

atire t fiows qUe wil1 not be usable for hi 9 h temper- 

Another approach^*®) for making measure- 
ments close to walls with the fringe LA is based 


on using the differences of the polarization prop- 
erties of the light scattered from particles and 
from solid surfaces. However, the success of this 
technique depends on the assumption of nonabsorbing 
spherical particles and on the details of the sur- 
face, such as roughness and the presence of 
deposits. 

In spite of its limitations, we feel that the 
fringe LA is suitable for many turbine engine com- 
ponent tests. For this reason, a detailed analysis 
of the fringe LA was conducted with the purpose of 
defining its limits of applicability. This paper 
addresses the problem of optimizing the optical 
system to achieve the best measurement of the 
velocity of a single seed particle. Note that 
this is only one part of the overall problem of 
the rapid, accurate determination of the flow 
parameters in turbine engine components. Other 
factors must also be considered, such as statis- 
tical biasing and proper data reduction methods. 

In this paper, the fringe-type laser anemom- 
eter is analyzed using the Cramer-Rao bound for 
the variance of the estimate of the Doppler freq- 
uency as a figure-of-merit. The noise from wall 
scatter is calculated using the wall bidirectional 
reflectivity and the irradiance of the incident 
beams. Mie scattering theory is used to calculate 
the Doppler signal . 

Similar work has been reported^, 10) for 
evaluating the signal-to-noise ratio for measure- 
ments close to walls. That work is extended here 
by treating the particle scattering using Mie 
scattering theory (instead of assuming a constant 
cross section) and by characterizing the wall 
reflectivity by the bidirectional reflectivity 
(instead of assuming a constant value for the 
reflectivity). Thus the analysis presented here 
can be used to evaluate the effect of various wall 
treatments and coatings which have reflectivities 
that are a function of both the incident light 
direction and the observation direction. 

In addition, a procedure is presented for 
determining the optimum aperture mask for the probe 
volume located a given distance from a wall. The 
use of circular aperture masks to decrease probe 
volume length has been recognized(^) and used in 
LA systems for compressor and turbine studies^ 3-6 ). 
However, the procedure presented here determines 
the optimum mask, which is not necessarily 
circular. 

Because the counter-type processor is usually 
used with fringe LA's for high velocity flows, the 
relation between the Cramer-Rao bound and the re- 
quired signal-to-noise ratio for counters is 
discussed. 

Finally, numerical examples for a coaxial 
backscatter LA are presented to illustrate the 
procedure for determining the optimum aperture 
mask. 


Theory 


In this section the Cramer-Rao lower bound 
for the variance of the estimate of the Doppler 
frequency will be obtained in terms of the param- 
eters of the Doppler burst. The parameters related 
to the signal will be determined using Mie scatter- 
ing theory. The noise parameter will be determined 
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by evaluating the laser light scattered by walls 
located near the probe volume. Finally, the per- 
formance of the counter-processor will be related 
to the Cramer-Rao bound. 

Doppler Burst Signal 

A particle passing through the probe volume 
of a dual-beam fringe- type laser anemometer 
(fig. 1) with the indicated trajectory (described 
by 5 , the distance of closest approach to the 
probe volume axis, and by $f, the angle be- 
tween the flow direction and the fringe normals) 
will generate a Doppler burst signal (assumed to 
be detected with a photomultiplier tube) 



the jo that maximizes p(nia) for a given set 
of measurements £.) 

For Poisson statistics, the probability of 
nj< counts in one time interval is 

-<i\> n. 

P(n k ) = e <n k > K / (5) 

where the expected value <n|<> is given by 

Eq. (1). The elements of Fisher information matrix 

can then be written^) 

ij - v<v 8a i 3a j ] 

If the parameters are uncorrelated, the vari- 
ance of the estimate of one of the parameters is 
given simply by the inverse of the corresponding 
diagonal element the Fisher information matrix; 
i.e., by the variance that would apply if only 
one parameter were unknown. 


where <n^> is the expected number of photo- 
electron counts in the time interval (tfc - At/2, 
tk + At/2), t„ is the time at the peak of the 
Gaussian envelope, fg is the "Doppler" frequency, 
fg is the frequency difference between the two 
incident beams, V is the visibility, N 0 is the 
constant background noise, Nf « 2 w Q /s is the 
number of fringes, (w 0 is the waist radius of the 
incident beams and s is the fringe spacing) and 


A = 


V 


-2(«/w q ) 


( 2 ) 


is the peak value of the pedestal, where A 0 is 
the peak value if the trajectory passes through 
the center of the probe volume. The trajectory is 
assumed to be near the focal plane and to have no 
significant component along the optical axis. The 
particle is assumed to have constant velocity as 
it passes through the probe volume. 

Cramer-Rao Lower Bound 


The unknown parameters in the Doppler burst 
can be denoted by the vector & = (fg, <j>f, t Q , A, 
il», N 0 , V). The measurement can be described by the 
vector £ = (ni, n?, nj, ..., njJ with elements 
that are the pnotoelectron counts for k equal time 
intervals At. 

The Cramer-Rao lower bound for the variance of 
the estimate of the parameter «j is'^J 

-ftty = (x -1 ),, (3) 


nbj > l/r 11 (7) 

Here, we assume that the parameters are only weakly 
correlated, so that this much simpler relation for 
the Cramer-Rao bound may be used. 

If we assume that the background noise is 
greater than the peak signal, and that the number 
of fringes is large compared to one, we can let 
At ♦ 0 and convert the sum in Eq. (6) to an 
integral, which can be evaluated to give a lower 
bound for the variance of the estimate of the 
Doppler frequency 



> 


64 

7T 5/2 


N f n J 
o D 


v 2 a 2 n 3 


(8) 


where we let N « Nf cos 6f; i.e, N is the number 
of cycles between the 1/e^ times of the envelope of 
the burst. Note that the lower bound is not de- 
pendent on the frequency shift fg. Equation (8) 
can be rewritten to give a lower bound for the 
relative error (corresponding to the Cramer-Rao 
bound) of measurements of the Doppler frequency 



where o f = [f (f D )]l/2 and where n 0 = N 0 /fp and 

a = A/fp are the expected noise and signal counts 
per Doppler period. 


where £ is the Fisher information matrix with 
elements given by 

2 

3 In p(n ia) 


The conditional probability function p(jn&) is 
the likelihood function, which expresses tTie prob- 
ability that given a set of parameters &, the set 
of observations jn will occur. (The maximum like- 
lihood estimate for the parameters a is given by 


A comparison (for a particular set of param- 
eters) between a numerical evaluation of (a. /f n ) 

t d u 

using the full Fisher information matrix (eq. (6)), 
and the analytical expression (eq. (9)), is shown 
in Fig. 2. This indicates that the assumptions 
used in the derivation of Eq. (9) are valid for 
n 0 /aV > 2. The simple expression for the error 
bound (eq. (9)) will be used in the remainder of 
this paper. 

The task now is — given a particular optical 
configuration, seed particle, test environment, 
and flow — to calculate the signal and noise. We 
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then can use the lower bound of the relative error 
given by Eq. (9) as a f igure-of-merit for this 
particular optical configuration. The optical 
configuration may then be optimized to find the 
maximum value of the lower bound (i.e., the small- 
est error in the measurement). 

Noise 

The background noise caused by laser light 
scattered from surfaces near the probe volume is a 
function of the optical configuration, the reflec- 
tance properties of the wall, and the distance of 
the wall from the probe volume. In addition to 
this direct wall scatter, there generally will 
also be a constant background noise flux caused by 
multiple reflections of the laser light, ambient 
light, and blackbody radiation from hot surfaces. 

The light flux through area dA at a point 
(x,y) on the aperture plane due to wall scatter 
and background radiance is obtained in Appendix A 
by integrating the radiance L of the wall over the 
solid angle subtended by the image of the pinhole 
at the point on the aperture (fig. 3). This light 
flux, expressed in terms of photoelectron counts 
per Doppler cycle is (eq. Al) 

n(x,y)dA *= 

[ L r (x ’ y; v* r ) + L b ] dn cos 6 dA 0°) 


For purposes of calculation, the aperture 
plane is assumed to be at the plane of the focusing 
lens instead of the position shown in Fig. 1 (the 
mask will have the same effect in either location). 
Also, for this calculation, we treat the laser 
light as if it were incoherent. This means that 
interference effects, such as speckle, are not 
taken into account. But the results obtained do 
represent the expected values. 

Signal 


polarized scattering is stronger than the "s" 
polarized scattering. Also, it can be seen that 
the phase angle variation is much larger for the 
larger particle. This phase variation is important 
because it leads to a reduction in the visibility 
of the Doppler burst when the signal is integrated 
over the aperture. 

Evaluation of Cramer-Rao Lower Bound 

The relations derived above for noise and 
signal are now used to evaluate the relative error 
given by Eq. (9) for a given optical configuration, 
seed particle, and test environment. 

By integrating over the clear aperture (not 
including light blocked by the aperture mask), we 
obtain the total number of counts for the two quad- 
rature signals, the pedestal, and the background 
noise: 

S ° ~ f s(x,y)dA 

q Q = f q(*,y)dA 

f . . (12) 

P 0 = J p(x,y)dA 

n ° = j n(x,y)dA 

These quantities are identified with the corre- 
sponding quantities in Eq. (1) for the Doppler 
burst signal ; i.e., 

A/f 0 = a = p 0 (13) 

is the pedestal component, 



is the component modulated at the Doppler 
frequency. 


For the geometry shown in Fig. 4, the signal 
is calculated in Appendix B using Mie scattering 
theory. The light flux through differential area 
element dA located in the aperture plane at posi- 
tion (x,y) consists of the following: the quadra- 
ture components of the Doppler modulated signal, 
s(x,y) and q(x,y); and the pedestal component, 
p(x,y). These are expressed in terms of the number 
of counts/Doppler period as 


s(x,y)dA = (2I/F 2 )(nX/hcf n ) cos 3 e CT O cos 6 dA 

u Tri a & <*B 

ct , p 

q(x,y)dA = (21 /F 2 ) ( n x/hcf n ) cos 3 e a e sin 6 dA 

u TTa a0 <*3 

a , p 

p(x,y)dA = (2I o /F 2 )(nx/hcfp) cos 3 e I sdc ( 11 ) 

The calculation takes into account both the 
variation in the amplitude and in the phase of the 
scattered fields as a function of scattering angle. 
Examples of the variation of the differential 
scattering cross section and phase angles are shown 
in Fig. 5 for 0.5, 1.0, and 1.5 ym diameter poly- 
styrene latex particles (refractive index = 1.59). 
The figure shows that for near backscatter the "p" 



is the visibil ity, and 

+ = tan' 1 (q 0 /s 0 ) (16) 


is the phase. 

For the background noise dominated case (for 
which Eq. (9) is valid), the Cramer-Rao lower bound 
for uncertainty in the measurement of the Doppler 
frequency is thus 



Counter-processor 

The above analysis gives a lower bound for the 
uncertainty in the measurement of the velocity of 
a single particle (eq. (17)). In practice, this 
error bound cannot actually be achieved. In par- 
ticular, the counter-processor exhibits a threshold 
at a signal-to-noise ratio (SNR) of about 10. C 14 ) 
For SNR > 10 the counter can be expected to give 
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good results, but below this level the measurement 
error increases rapidly. As in Ref. 14, we define 
the SNR as the ratio of the mean signal power at 
the peak of the Doppler burst to the noise power 
due to wall scatter and background 


But, because of symmetry, only one quadrant need 
be considered with 100 area elements. 

The error bound can be expressed in terms of 
these 100 elements as 


SNR (AV ) 2 S o * q o 

SNR = -ZTil^TT - 4n 0 (Af/f 0 ) 


(18) 



where Af is the noise bandwidth. (The noise 
bandwidth is determined by the low-pass and high- 
pass filters in the processor.) 


Comparing this expression (eq. (18)) for the 
SNR with Eq. (17) allows us to express the SNR in 
terms of the relative error (o^ /f D ) given by 

Eq. (9). 0 



(19) 


Even with SNR > 10, the counter will not pro- 
vide measurements as accurate as indicated possible 
by the Cramer-Rao bound because of other error 
sources such as the finite clock frequency and the 
effectiveness of validation circuitry. 


where 


n k = n ^ x k’ y k^ iA * s k = ( x k ,y k^ aA » n k = n ^ x k ,y k^ aA 

(21) 


The noise elements nu and signal elements, s^ 
and q^, are evaluated by numerically integrating 
Eqs. (12). Standard Mie scattering algorithms are 
used for the differential cross sections. 


To determine the optimum aperture mask a 
numerical procedure is used to eliminate the set 
of areas that will maximize the quantity 


Q = 


(ls k ) 2 + Qq k ) 2 

£"k 


( 22 ) 


Numerical Examples 

Two examples of calculating the Cramer-Rao 
bound and selecting the optimum aperture mask are 
presented in this section. The following pa- 
rameters are used for both examples: the fluid 
velocity is 500 m/sec; the seed particles are poly- 
styrene spheres with 1 pm diameter and real refrac- 
tive index equal 1.59; no Bragg shift is used; the 
flow is perpendicular to the fringes; the laser 
power equals 0.1 W at 0.5145 pm wavelength; the 
beams are linearly polarized with the electric 
vector normal to the plane formed by the beams 
(along the y axis); and system efficiency including 
losses and PMT quantum efficiency is 10 percent. 

The wall is normal to the optical axis, is located 
1 mm from the center of the probe volume, and is a 
Lambertian reflector with reflectivity equal 1/*. 

For the first example, which represents a low 
performance optical configuration: the focusing 
lens has a 200 mm focal length and a 50 mm diameter 
(f/4); the pinhole image diameter is set equal to 
125 pm; the number of fringes is 12.5; the fringe 
spacing is 10 pp; and the uniform background radi- 
ance is 3. 1x10" 5 W/mm 2 /sr. 

For the second example, which represents a 
high performance optical configuration: the focus- 
ing lens has a 250 mm focal length and a 100 mm 
diameter (f/ 2.5); the pinhole image diameter is 
set equal to 70 pm; the number of fringes is 10; 
the fringe spacing is 5 pm: and the uniform back- 
ground radiance is 3.8xl0" 3 W/mm 2 /sr. Note 
that the background radiance is set higher for 
this example than it was for the first example. 

This was done to insure that the assumption that 
the noise is larger than the signal remains 
valid. Of course, this indicates a more severe 
environment. 

For these numerical calculations, the the 
aperture plane is divided in 400 equal areas. 


The procedure used is to calculate Q using 
all the area elements, then eliminate each of the 
areas, one at a time, and recalculate Q. The 
area element whose elimination gives the greatest 
increase in Q is included in the aperture mask. 
This procedure is repeated until no further reduc- 
tion in the error is possible. The resulting mask 
is the optimum for the particular distance from 
the probe volume to the wall. 

The calculated masks for the two examples are 
shown in Fig. 6. The error bound (o^ /f^) ■ 0.05 

and 0.002 for the f/4 and f/2.5 examples, respec- 
tively. (The noise bandwidth Af was set equal 
to the Doppler frequency fp.) The corresponding 
SNR = -7.1 dB and 25.8 dB, which means that a 
counter processor should be a suitable processor 
for the f/2.5 example, but not for the f/4 example. 
The mask covers 57 percent of the lens aperture 
for the f/4 optics and 51 percent for the f/2.5 
optics. Also shown on Fig. 6 is the order that the 
grid elements are blocked in the optimization 
procedure. 

Figure 7 illustrates how the SNR improves as 
the masks are built up. Note that the SNR improves 
less quickly as the masks approach their optimum 
configurations. This is particularly noticeable 
for the f/2.5 optics where the SNR would only be 
slightly decreased by using a mask that covers 20 
percent of the lens aperture instead of the optimum 
mask which covers 51 percent. 

Figure 8 shows the SNR as a function of dis- 
tance from the probe volume to the wall (the SNR 
was calculated for the optimum mask at each loca- 
tion). Note that the SNR reaches a constant value 
at large distances because of the inclusion of a 
constant background radiance. This figure illus- 
trates the improvement of the f/2.5 optics compared 
to the f/4 optics. With f/4 optics, measurements 
can be made to about 2 mm from the wall; while, 
with f/2.5 optics, measurements can be made to 
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about 0.5 mm. Of course, these results are depen- 
dent on the parameters chosen for the examples. 

In particular, the reflectivity of the wall has a 
large effect on the SNR for measurements near the 
wall. 

Discussion 

The optical design of a fringe-type laser 
anemometer for a particular experiment requires 
several steps to take into account all contraints 
imposed by the experiment and by the LA system. 

The initial step is to select a seed material that 
will follow the flow accelerations and will survive 
the flow environment. It is usually desirable to 
use monodisperse particles; however, this may be 
difficult to achieve in practice. 

The next step is to select the fringe spacing 
and number of fringes. The fringe spacing must be 
such that the maximum velocity will not result in 
Doppler frequencies greater than the maximum fre- 
quency response of the photodetector and the signal 
processor. Furthermore, the minimum fringe spac- 
ing may be limited by the available solid viewing 
angle. (The fringe spacing is inversely propor- 
tional to the beam crossing angle and hence to the 
beam separation at the focusing lens.) 

Another factor that can affect the selection 
of the fringe spacing is the phase cancellation 
that occurs when the scattered light is integrated 
over a large aperture. This can be a significant 
problem when the fringe spacing is on the order of 
the particle size. In fact, this phase cancella- 
tion is the basis for particle size measurements 
based ^tj^ visibility of the Doppler 

The selection of the number of fringes is 
based on several considerations. For counter pro- 
cessors with a given time interval measurement 
accuracy (typically, 1 ns.), the number of cycles 
used must be large enough to result in an accept- 
able error on a single particle measurement. Note 
that the acceptable single particle measurement 
error will be larger for highly turbulent 
flows. U/J 

The problem of statistical biasing, which 
occurs in turbulent flow, may also influence the 
selection of the number of fringes. For example, 
the biasing caused by flow angle fluctuations^^/ 
is reduced by using a larger number of fringes. 

Also, the use of frequency shifting can reduce 
some biasing errors (but at the expense of increas- 
ing the required frequency response of the photo 
detector and signal processor). 

After the number and spacing of the fringes 
have been established, the field stop (pinhole) 
can be sized to match the probe volume diameter. 

The focusing lens is selected based on the required 
working distance and available solid angle. It is 
generally advantageous to select as fast a lens as 
is practical. 

With the optical system specified, the pro- 
cedure presented in this paper for determining the 
optimum aperture mask can be applied. The mask 
should be optimized for the closest anticipated 
measurement position to the wall. Based on the 
calculated SNR, the suitability of a counter pro- 
cessor can be verified. 


Note that to apply this procedure the bidirec- 
tional reflectivity of the surfaces in the test 
region must be known. Of course, the surfaces 
should be treated to minimize the reflectivity. 
(This can be accomplished either by absorbing the 
light or by reflecting the light in a "safe" direc- 
tion.) For measurements near the window, the re- 
flectivity of the window must also be known. These 
reflectivities will generally have to be obtained 
experimentally. It should be realized that the 
reflectivities can change during a test run as 
contaminants accumulate on the walls and windows. 

In addition, the background radiance due to flames, 
blackbody radiation, and multiple reflections of 
the laser light must be estimated or measured. 

Finally, it should be noted that the assump- 
tion of spherical seed particles that was used for 
the Mie scattering calculations may introduce some 
error if the particles are not, in fact, spheres. 
And, if polydisperse particle size distributions 
are used, the analysis should be modified to take 
this into account. 

Concluding Remarks 

Using the analysis presented in this paper, 
the fringe-type laser anemometer can be evaluated 
in terms of its best possible performance based on 
the Cramer-Rao bound for the estimate of the vari- 
ance of the measured Doppler frequency. The rela- 
tion between this lower bound and the performance 
of counter-type processors points out the potential 
performance gain that should be achievable using 
more optimal processors for low signal-to-noi se 
ratios. 

APPENDIX A - CALCULATION OF NOISE 

The irradiance at a point (x,y) on the aper- 
ture plane due to wall scatter is obtained by in- 
tegrating the radiance L of the surface over the 
solid angle subtended by the image of the pinhole 
at the point on the aperture (fig. 3). For pur- 
poses of calculation, the aperture plane is assumed 
to be at the plane of the focusing lens instead of 
the position shown in Fig. 1. The light flux 
through area dA at point (x,y) on the aperture 
expressed in photoelectron counts per Doppler 
cycle is 

n(x,y)dA = 


(hc f D ) f [ L r ( x >y; ® r >* r ) + L b Jdn cos e dA ( Al) 

The radiance due to single scattering of the inci- 
dent beams from the wall is 

L r (x,y;e r> $ r ) = p , (x',y , ; 0 i ,*.;e r , 4 , r ) I.fx'.y') 

(A2) 

where I(x',y') is the irradiance of the beams at 
the surface and p' is the bidirectional reflec- 
tivity. Note that p' is a function of both the 
direction of incident and the reflected light. 

If the wall is normal to the optical axis, 
the irradiance of the Gaussian beams at the surface 


6 





l ~ 2(xf + y 2 ) 




W^Zj) 


- 2<x 2 + y 2 ) 


w 2 (Zg) w 2 (Zg) 


(A3) 


where (x‘, y 1 ) is the point on the wail along the 
ray through point (x,y) in the direction (e,*). 
Also, (xj, yj, zx) and (x2, y2, z z) are the coordi- 
nates of the point on the wall in coordinate system 
with zx and Z2 along the beam propagation di- 
rections and are given by 


X l,2 


COS B 

0 

+ sin b 


y l,2 

= 

0 

1 

0 


z !,2_ 


+ sin b 

0 

COS B 



with d being the distance from the focal point 
to the wall and w(z) the beam radius given by 



APPENDIX B - CALCULATION OF SIGNAL 


The signal is calculated using Mie scattering 
theory, (19) w hich takes into account both the 
variation in the amplitude and phase of the 
scattered fields as a function of scattering angle. 
The incident beams are located in the x-z plane 
and have a crossing angle 2$ (fig. 4). 

The electric fields of the incident beams 
near the waists are represented as plane waves 


1,2 " e l ,2 l 01,02 


T2,f lj2 t-ik 1>2 .r 


(Bl) 


where the field amplitudes are Eqi q 2 = (2Z p ) 1 I 0 ; 
with I 0 being the peak irradi ance’of each 6eam at 
its waist, which is related to the laser power and 
waist diameter by I 0 = P^/ttw^. The electric 
fields of the beams are assumed to be plane polar- 
ized with unit vectors ^x and ej?. 

We define the following unit vectors: 


2 = (k s x k lj2 ]/|k s x ki f2 | 

(B2) 

Pi, 2 = < k l,2 x $ l,2 )/k 

(B3 ) 

q i>2 = (fc s x ^i,2 )/k 

(B4) 


where ^x $2 a X e normal To the scatter- 

ing planes; p^ and P 2 are normal to the 
propagation directions of the incident beams and in 
the scattering planes; and, qx and q? are 
normal to the scattered wave propagation direction 
and in the scattering planes. 


The incident fields can be expressed in terms 
of the components in the scattering plane and nor- 
mal to the scattering plane 

g l,2 = E 1s,2s *1,2 + E lp,2p ^1,2 (B5) 

The scattered fields are given by Mie scatter- 
ing theory for spheres (*9) 

-1, S 2 = 1 6 krt [ E 1 s,2s S s ^®1,2^ ^1,2 

+ E lp,2p S p (0 1,2^ ^1,2^ (B6) 

where R is the distance from the scatterer to the 
observation point, and S § (e) and S p (e) are the 
usual Mie scattering amplitude functions with the 
scattering angles given by 

cos e lj2 = • k s )/k 2 ( B7 ) 


The irradiance of the scattered light at the ob- 
servation point is 



This is composed of the sum of a slowly varying dc 
component, I s dc* corresponding to the pedestal 
of the Doppler burst and an ac component, I sac , 
modulated at the Doppler frequency. These can be 
written as 


- >. [‘W « s <«i> * « P <»i> 

+ (e 2 *s 2 ) 2 (e 2 ) + (e 2 -p 2 ) 2 o p (e 2 )j 

= l o 


(B9) 


a=S,p 

B=S,p 

(cos 6 cos 2irf n t - sin 6 n sin 2irf n t) (BIO) 

aB D aB D ‘ v ' 


Note that I $ac is expressed in a quadrature 
representation form. Because the phase of I sac 
varies across the aperture, the detected signal 
will, in general, be less than the signal that 
would be detected if the phase were constant. 

The Mie complex amplitude functions can be 
written as 


i 6 (e) 

S„(e) = |S o (e)| e “ , a = s.p (Bll) 

The differential cross sections for "s" and "p" 
scattering are 

o a (e) = U/2n) iS a (e)i 2 , o = s.p (B12) 

In Eq. (BIO) we also introduced the following 

°ss = [o s (e l> °s (e 2^ 1/2 <V*1> (^2^ 2 ) <V $ 2> 
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°PP = [a p (6 l) °p( 9 2^ 1/2 44> 44) (^4) 

°s P - [o s (e i) y e 2>] 1/2 44 ) 44 ) 44 > 
v = o P 4 > ° s ( e 2)) 1/2 44 ) 44 ) 44 ) 

*a B = S a (9 l> - S b { 9 2 ) * a ’ 8 = S *P (B13) 
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/-PINHOLE IMAGE 



Figure 1. - Optical configuration of coaxial backscatter fringe- 
type laser anemometer. 



N 0 /aV 


Figure 2. - Ratio of approximate value (eq. 9) to exact 
value (eq. 3) of lower bound for relative error of Doppler 
frequency versus ratio of background amplitude to sig- 
nal amplitude; Parameters used in calculation are: 
number fringes, N = 10; signal amplitude, aV = 50 
counts/ Doppler cycle? visibility, V - 1 ; and, phase, 

4= 0. 



Figure 3. - Geometry used for evaluation of wall scatter noise. 



Figure 4. - Geometry used for evaluation of Mie scattering. 










(a) f/ 4 Optics. 

Figure 6. - Optimum masks for probe volume located 1 mm from wall Grid 
elements are numbered in the order that the mask is built up. 


y 

♦ 



(b) f/2. 5 Optics. 
Figure 6. -Concluded. 



SNR, dB 



(b) f/ 2. 5 Optics. 

Figure 7. - Signal-to-noise ratio (SNR) versus accumulation 
of mask area elements. 
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(a) f/ 4 Optics. 


Figure 8. - Signal-to-noise ratio (SNR) 
for optimum mask versus distance of 
probe volume from wall. 



(b) f/2. 5 Optics. 


Figure 8. - Concluded 
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